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Abstract 

A hypothesis is put forward regarding the function ^(x) which describes the dis- 
tribution of twin primes in the set of natural numbers. The function TT2(x) is tested by 
evaluation and an empirical vr^x) is arrived at, which is shown to be highly accurate. 
Several other questions are also addressed. 
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1 Introduction 

In 1923, Hardy and Littlewood proposed a hypothesis for the distribution of twin primes on 
the interval (l,x) pQ: 

7T 2 (x) ~ 2 TT ( 1 - 7 ^ n • (1-1) 

1 ; f = \\ (p-i)V 

Later on, the following expression was put forward [2J: 

vr 2 (x) ~ 2 FT ( 1 - — L — ) -^— 2 T\ (1.2) 
i = \\ (P-1)V (\nx) 2 ^\p-2 



The asymptotic representations (11 .11) and (11.21) of the very important function ^(x) are too 
complicated to be used in practice. In this article I propose a new law for the distribution of 
twin primes among the naturals in the form of a much simpler 7^(2), based on composition 
with the function ir(x). 

2 Results 

Recall that tt(x) [^(x)} is the number of primes (twin primes) not larger than x. The 
following hypothesis is proposed for the distribution of twin primes in the set of all naturals: 
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Hypothesis 2.1. Twin primes are distributed among prime numbers in the same way that 
primes are distributed among naturals. In other words, 

TT 2 (x) = IT (tt(x)) . (2.3) 

Table [□ gives some values of the functions 7r 2 (a;) and tt (tt(x)) for x < 10 6 . The values of 
7T2(x) were computed according to Lehmer's tables [3]. From the results in Tabled] it is safe 
to say that the ratio of 7r 2 (a;) to n (7r(a;)) is either exactly one or differs from unity by some 
negligibly small amount. 
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Table 1: Testing hypothesis (12. 3 1) for a few selected values of x < 10 6 . 
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Next, we will require the upper and lower bounds for 7r(x) given in [I]: 



< 7r(ar) < — — . (2.4) 



3 In x 5 In a; 

Theorem 2.1. For all x > 5 for which ft2.4\ ) holds, we have 

A < vr 2 (x) = tt{tt{x)) < B, (2.5) 

where 

A - iX 



91nx [In a; — ln(lnx) — In 1.5] ' 

64x 

25 In x [In x — ln(ln x) + In 1.6] 



Proof. The function f(x) = x/lnx is monotonically increasing for x > 3. From (12. 4p we 
have 

27r(x) . . , . 8tt(x) 

- <7T 2 (x) = 7T 7T X < , r / N1 ■ 2.6 

31h[tt(x)] K ' K K " 51n[7r(x)] v ' 

If we now consider the expression as a function of 7r(x), we can see it is also mono- 

tonically increasing. In our case, ir(x) > 3 as x > 5 (by theorem requirements). Taking this 
into consideration and using the right-hand side of inequality (I2.4p 

/ \ , , ,s 8vr(x) 8^ 64x 

n 2 (x) = Ti Mx)) < ; r / \ n < % o \ = — : f; r-^ ; ; r = 

w v K " 51n[7r(x)] 51n(g^) 25 lnx [lnx - ln(lnx) + In 1.6] 

Similarly, by using the left-hand side of inequality (12. 4 \ we obtain the lower bound: 

V; V V ;; 31n[vr(x)] 3^(3^) 91nx [lnx - ln(lnx) - In 1.5] 
Thus, from ( 12. 6ft we have 

. 27r(x) , . . , 8n(x) 

A < sbFM < * 2W = " < FtaFwI * B ' 

precisely inequality ( 12. 5ft . as required. □ 
In Table |2] we check inequality (12. 5 ft for several values of x. 

Next let us look at the density of twin primes among the primes. 
Theorem 2.2. Almost all primes are not twins, so 

tt 2 (x) = o(tt(x)). (2.7) 
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Table 2: Testing inequality ( 12. 5 j) for a few selected values of x < 10 6 . 



Proof. Assume that hypothesis ( 12. 31) is true. Then, denoting y = tt(x), we have 

7T 2 (x) 71 Mx)) , . . 

< -^f = = 7T (y )y. 

We can find an upper bound for n(y)/y by the sieve method, taking the set {y} to contain 
no repeated values. 

Let ip{y,r) be the number of naturals no larger than y and not divisible by any of the 
first r primes Pi, P2, . . . , P r . Then 

<p{y,r)= Yl Kd)[y/d\, (2.8) 

d\P 1 P 2 ...P^ V y ) 

where fi(d) is the Mobius function and d\P\ . . . P r means all d not divisible by Pi to P r . It 
is clear that 

7r{y)<ip{y,r)+r. (2.9) 

We next drop the floor operator in ( 12.81) . and note that there are 2 r terms being summed. 
This means that the resulting expression has an error no larger than 2 r , and by ( 12. 9ft we 
subsequently get 

Ay) < E Kd)[y/d\+r< y x ^r + r + r 

d\PlP2...P^) d|PlP2...P w(V 5) 

= y J] (l-P^) + r + 2 r <y J] (l-p- 1 )+2 ? ' +1 , 

P<P r P<P r 
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because r < P r < 2 r . Furthermore, from the inequality 

nti-p-y^in 



we find 



P<x 



<v) <r^ + i r+1 < r~ + 2 r+1 . 

in Fr in r 



Choose r = clny, c < l/ln2. Then 2 r < y and 



^ in [c in y] ^ in c + in [in y] ^ 



where cln2 < 1. Dividing through by y, we get 



n(y) 1 2 
< — < : r + 



y lnc + ln(ln?/) ' y 1 cln2 

As y — > oo, the right-hand side of the above goes to zero, which implies the validity of 
equation (12. 7p . □ 

Corollary 2.1. Since we know that tt(x) = o(x), from Theorem \2.S\ it follows that ^(x) = 
o(x) . In fact, 

nm ^ = ^ ^x) Um vr(x) = Q 

x-»oo a; a;-»oo 7I"(x) x->oo X 

Next, we move on to construct an empirical function for the law of distribution of twin primes. 

Denote by rjp the density of primes in the reals, and by t]pp the density of twin primes 
in the primes, i.e. rjp = ir(x)/x and r\pp = 7r 2 (x)/7r(x). Based on 7r(x) = o(x) and (12.71) . the 
densities rjp and r/pp go to zero as x — > oo, but the ratio 

h=^ (2.10) 
Vp 

remains bounded in a well-defined, constant interval (see Table [3]). We can obtain a rough 
estimate of an upper bound for h > 0; for this we need the inequality tc(x) > x/ lnx and the 
right-hand side of (12. 5ft . We get 

X7r 2 (x) 64x 2 



[tt(x)] (^L) 251nx(lnx-ln(lnx) + In 1.6) 

2.56 lnx 2.56 lnx _ _ „ 
< ; ^—r, r < 5.12. 



lnx — ln(lnx) + In 1.6 lnx — ln(lnx) 

Thus, < h < 5.12. This fact allows one to construct an empirical function ^(x) for the 
number of twin primes on (2,x). As is evident from Table [3], ^(x) defined below is rather 
accurate. 
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We obtain ir%(x) / ir(x) = h c ■ n(x)/x, leading to 



x 

where h c = 1.325067 ... - the mean value of h for x < 10 6 and we round the right-hand side 
of (12. lip to the nearest integer. 

In Table [3] I test the accuracy of 7r|(:r) for 50 < x < 10 6 . Nevertheless, (12. lip is appli- 
cable for x > 10 6 , too. For example, there are 183, 728 twin primes less than or equal to 
x = 37 • 10 6 , while vr^x) = 183,463 which gives a relative error of 6 = 0.0014 (see Table [3]). 
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Table 3: Testing expression (12. lip for a few selected values of x < 10 6 . 
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